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Abstract 

We analysed the ground state energy of some dimerized spin-i transverse XX and Heisenberg chains 
with Dzyaloshinskii-Moriya (DM) interaction to study the influence of the latter interaction on the spin- 
Peierls instability. We found that DM interaction may act either in favour of the dimerization or against it. 
The actual result depends on the dependence of DM interaction on the distortion amplitude in comparison 
with such dependence for the isotropic exchange interaction. 
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The spin-Peierls instability is known as a magnetic analogue of the conventional Peierls instability in electron- 
phonon systems. A uniform quantum spin chain at low temperatures may become unstable towards dimerization 
owing to the interaction with lattice degrees of freedom. This occurs because the dimerized lattice distortion 
lowers the magnetic energy by a greater amount than the increase in the elastic energy due to deformation. 
Starting in the 70s with organic compounds exhibiting spin-Peierls transition, the interest in the spin-Peierls 
instability of quantum spin chains was renewed with the discovery of the inorganic spin-Peierls compound 
CuGe0 3 in 1993 |, g]. To model appropriately the spin degrees of freedom of the spin-Peierls compounds 
the pure Heisenberg chain as well as its modifications, which include frustration or interchain interaction, are 
considered. As a rule, since those models represent quantum many body systems, only approximate results can 
be obtained. However, some generic features of the spin-Peierls systems can be illustrated in a simplified but 
exactly solvable quantum spin model, namely, the transverse XX chain || |], ^ |[ . 

In the present paper we discuss the influence of the Dzyaloshinskii-Moriya (DM) interaction on the spin- 
Peierls dimerization in the adiabatic limit. The presence of DM interaction for CuGe03 was proposed in Refs. 
^ [I [l0| in order to explain the EPR and ESR experimental data. The structure of DM interaction in the 
cuprates was examined in Refs. JTT| , PL The influence of this interaction on the ground state properties of 
the one-dimensional and two-dimensional Heisenberg models was studied in Refs. Besides, the 

multisublattice transverse XX chain with DM interaction was introduced in Ref. Jl^j , however, the spin-Peierls 
instability was discussed only in one limiting case (see below) in the absence of an external field. DM interaction 
was found to be present in a number of quasi-one-dimensional magnets (see, e.g., (Tt], [l^] and also [|l9) in which 
RbCoCl3-2H2 is described as a pure DM chain) and a study of one-dimensional DM Hamiltonians seems to 
be of great importance. 

Several mechanisms which may destroy the dimerized phase in the Heisenberg chain were discussed in the 
literature, in particular, an external field or an Ising anisotropy. XX anisotropy also suppresses the dimerized 
phase although does not destroy it completely. Therefore, initially one may expect that the appearance of 
an anisotropy in form of the DM interaction could act against dimerization. Although we shall find that the 
Heisenberg chain with increasing of DM interaction may become similar to the XX chain both enhancing or 
suppressing of the dimerized phase are possible depending on the details of the distortion dependence of DM 
interaction. 

In our study we follow the idea of Ref. || and compare the total ground state energy of the dimerized and 
uniform chains in the presence of DM interaction. First we demonstrate that in some cases the DM interaction 
can be eliminated by a spin coordinate transformation resulting in a model with an anisotropic exchange 
interaction. This observation permits us to study rigorously the influence of DM interaction on the spin-Peierls 
instability in the transverse XX chain using the exact results for thermodynamic quantities of the regularly 



alternating transverse XX chain obtained recently with the help of continued fractions p0[ . (The approach 
exploiting continued fractions in contrast to the approaches used in previous works ]3|, ^, ||, ^) allows one to 
consider in a similar way not only the dimerized lattice but also more complicated lattice distortions.) Further, 
we discuss the case of the Heisenberg chain with DM interaction using exact diagonalization of finite chains. The 
exact analytical findings for the XX chain are helpful for the interpretation of the finite-chain results for more 
realistic Heisenberg chain. In the present study we are interested in describing the generic features originated 
by DM interaction and therefore no compound specific parameters are considered. 

To begin with, we consider a nonuniform chain of N — > oo spins | governed by the isotropic Heisenberg 
Hamiltonian with DM interaction 

H = + ^2 («/ n (s n • s„+i) + D„ ■ [s„ x s„+i]) . (1) 

n n 

Here VL n is an external field at site n, and J n and D„ are the isotropic exchange interaction and the antisymmetric 
anisotropic exchange interaction or DM interaction between the neighbouring sites n and n + 1, respectively. 

Extending the spin coordinate transformation used in Ref. |2l| for models with nonuniform J n and D ra first 
we show how the terms with cross product can be eliminated from Hamiltonian (Q) in special cases. If D„ 
has only one nonzero component D z n one may perform local rotations about the z axis introducing new spin 
operators = s* cos (j) n + s v n sin </>„, s v n ' = —s* sin (f> n + s v n cos </>„, <b n — ipi + . . . + <p n -i, tan ip m — in terms 
of which iMi becomes |23l 



H = J2 *W + E f V J « + D n («+x + + ■ (2) 

Note, that such transformation can be also applied to the transverse XX chain in which z spin components do 
not interact. Assume further that D„ has only one nonzero component D^. Then the local rotations must be 
performed about the x axis, i.e. = s^cos0 n + s^sin0„, s z n — — s^sin0„ + s^cos</>„ with tan^? m = -j 22 - 
resulting in 



H = ^n (Shl«' + COS«') + J] + \J ' <% + D? + 



(3) 



Acting similarly for D„ = (0, £^,0) one finds that the transformed Hamiltonian is given by (^) with the 
replacement — > D^, s xl — ► s v ', s v ' — > — s x> . In the case when D„ has more than one nonzero com- 
ponents the described elimination can be performed if the orientation of D„ (but not necessarily its value 



D n = y D^ 2 + D % n + D^ 2 ) is site independent, i.e. ^f - do not depend on n. (Obviously, the partic- 
ular site independent orientations of D„ reproduce the cases discussed above.) Really, in such a case we start 
from the global transformation of coordinate system with the Eulerian angles 0, 9, ip putting tan^ = -pg-, 9 = ^, 
tan ib — . " P " and then perform the mentioned above local rotations about the x axis with tan 09„, = ^ m - 



finding as a result 



H = E fi » ( + " At ( cos «' - sm M') 



E + V J n + A, 2 + <'<+/) ■ (4) 



It should be noted that the described transformations were applied to open chains. For cyclic chains they yield 
the presented expressions (||) - (Q) at least up to the boundary term. Evidently, the thermodynamic properties 
of the initial and transformed Hamiltonians are identical. Due to these transformations (Eqs. (|J) - (|J)) one 
is able to exploit the broad knowledge on anisotropic Heisenberg chains and XX chains. In what follow we 
shall use the thermodynamic equivalence of the initial and transformed Hamiltonians both in the analytical 
treatment and numerical computations. 

We proceed considering the transverse XX chain with DM interaction having only z component 

H = E « + E J « + S n S n+l) + E D n " «+l) ■ (5) 

n n n 

For the Hamiltonian (||) we are able to perform rigorous analytical calculations since after the Jordan- Wigner 
transformation it reduces to noninteracting spinless fermions. As it has been already mentioned thermodynamics 
of the model given by (||) is the same as of the transverse XX chain (without DM interaction) with exchange 



interaction y J 2 + D^ 2 . The thermodynamic properties of the latter model with regular alternation in bonds 
and fields having finite period p have been examined recently with the help of continued fractions p0[ . To 
study the spin-Peierls dimerization in the adiabatic limit we need the ground state energy of a spin chain with 
period p — 2, i.e. with the sequence of parameters Sli J\D\ f^^-Df^i^i-^i^^-Df • ■ ■ • Moreover, we assume 
the following reasonable dimerization ansatz J\ = J(l + S), J2 = J(l — S), D± — D z (l + k5), D\ = D z (l — kS), 
where < i5 < 1 is the dimerization parameter. It is argued that the directions of D-vectors are not changed by 
the dimerization (h], 23 . From Ref. @] we know that the dependence on the intersite distance of the isotropic 



exchange interaction and DM interaction may be different. This effect is described by the parameter k. Putting 
k = one has a chain in which D z does not depend on the lattice distortion, whereas for k = 1 the dependence 
of D z on the lattice distortion is as that for the isotropic exchange interaction J. The latter case with 57„ = 
was considered in ]l6|] . Besides, we bear in mind that, as a rule, the value of DM interaction is significantly 
smaller than the value of isotropic exchange interaction H . 

Further, we consider the case of zero temperature and look for the total energy per site £(S) which consists 
of the magnetic part eo(S) and the elastic part aS 2 , a > 0. From Ref. [^0| we know the exact expression for 
the magnetic ground state energy 



■■^ = ~w(* a 3 a )-5 |ni + n " l (5-?l- 



bi,2 = ^(f>i-^ 2 ) 2 + {h±h)\ 



h,2 = V J2 (! ± (5 ) 2 + Dz2 ( l ± M ) 2 > 
if bi < ||0 X +« 2 ], 



arcsinj ^;^'^ , if b 2 < + fi 2 | < bi, 



if i|fii + n 2 | < b 2 , 



where E(^,a 2 ) = <WT — a 2 sin 2 (/> is the elliptic integral of the second kind. We also seek for a nonzero 
solution 8* 7^ of the equation -4^ = that can be easily derived from (^). In what follows we consider the 
case of a uniform transverse field fii = il 2 = f^o > 0. In the limit J « 1 valid for hard lattices (having large 
values of a and corresponding to the experimental situation) one finds bi = I,b 2 = IKS with I = VJ 2 +D z2 
and N = ^jrir§^ ■ Instead of Eq. (|^) one then has 

e (5) = --E^,l-^5 2 )-n (l-t\ 7 (7 ) 

7T \ Z TT I 



if J < ft , 



/2_Q2 

arcsin^/ ja^,^) , if •/«<) < S2 < /, 

if n < ms, 



4> 



whereas the equation for 8* reads 

where F(ip, a 2 ) = \f , d ^ is the elliptic integral of the first kind. 

u Y 1 — a 2 sin 2 

Consider at first the case Oo = 0. After rescaling / — > J, ^ — > a, K<5* — > (5* one finds that Eq. (JsJ) is exactly 
the same as considered in Ref. || and thus 6* ~ 4 exp (— tw)- Thus for fe = 1 (N = 1) nonzero D z leads to 
an increasing of the dimerization parameter 8*, whereas for k = (H < 1) nonzero -D z leads to a decreasing of 
<5*. Let us pass to the case < Qq < I. Varying 8* in the r.h.s. of Eq. (@) from to 1 one calculates a lattice 
parameter a for which the taken value of 8* realizes an extremum of £(8) (^). One immediately observes that 
for j$ < 8* the dependence a versus 5* remains as that in the absence of the field, whereas for < 8* < 
the calculated quantity a starts to decrease. From this one concludes that the field ^j- — exp (— fror) makes 
the dimerization unstable against the uniform phase. The latter relation tells us that nonzero D z increases the 
value of that field for k — 1 and decreases it for k — 0. It is known Q that the increasing of the external field 
leads to a transition from the dimerized phase to the incommensurate phase rather than to the uniform phase. 
However, the former phase cannot appear within the frames of the adopted ansatz for the lattice distortions 
8i8 2 8 1 S 2 . . . , 8i + S 2 = 0. 

After the discussion of the limit 8 <C 1 we now present the results for arbitrary < 8 < 1 based on (^). In 
Figs. 1, 2 we plot the changes of the total energy £ (8) — £(0) (0) vs 8 and the nonzero solution 8* of equation 



q S ' = vs a, respectively, for various strengths of DM interaction. These results confirm that for k = 1 
DM interaction D z acts in favour of dimerization, whereas for k = against it. To understand the validity 
of the data obtained by exact diagonalization of finite chains that will be used below for the Heisenberg chain 
we present in Fig. 1 also the numerical results £{8) — £(0) vs 8 for N = 24 spins (open boundary conditions). 
Note, that although the finite chain results for N = 24 still overestimate noticeably the value of 8* and the 
depth £(8*) — £(0), however, and this is most important, they reproduce qualitatively correctly the influence 
of D z in both cases k = 1 and k = 0. The possible influences of D z are reproduced correctly even for shorter 
chains of N = 16, 20 spins. Moreover, with increasing the chain length from TV = 16 to 24 the numerical data 
approach the analytical ones valid for N — > oo. This expected tendency can be also seen by comparison curves 
1 for N = 24, N = 28 (the computation in this case becomes already very time consuming) and N — > oo in Fig. 
1. 

Let us turn to the Heisenberg chain with DM interaction. Assume that the vectors D n have the same 
orientation at all sites, e.g. in the z direction. (The assumption D„ = (0, 0, D^) does not lead to a loss 
of generality if Oo = 0.) In such case the influence of DM interaction on thermodynamics follows from a 
study of thermodynamic properties of the Heisenberg chain with anisotropic exchange interaction (Q). We 
immediately find the appearance of XX anisotropy, since the interaction between x and y spin components 
becomes \J J 2 {1 ± 8) 2 + D 2 (l ± kS) 2 (» \J J 2 + D 2 (l ± #8) for i5<1) whereas between z spin components it 
remains J(l±8). Restricting ourselves to small 8 we perform the redefinitions \J J 2 + D 2 = J', = 8', ^ = a'. 
As a result we come to the anisotropic alternating Heisenberg chain with xx and yy interactions J'(l ± 8') and 
zz interaction J'A(1 ± ff), A = -j=^=== and the expression for the elastic energy per site a'8' 2 . Any increase of 
D results in increasing of intersite coupling J' and increasing of XX anisotropy manifested by going of A from 
1 {D = 0) to (D — > oo). Hence, the Heisenberg chain should start to exhibit a behaviour inherent in the XX 
chain while D becomes large. Omitting the role of a change in the zz interaction, which becomes less important 
while D increases, and bearing in mind the corresponding analysis for the XX chain based on (^), we may 
expect that for £!o = the dimerized phase will be enhanced for k — 1 (since a' = a) and will be suppressed 
for k = (since a' > a). The numerical results for £(8) — £(0) vs 8 for the Heisenberg chain of N = 24 sites 
(open boundary conditions) shown in Fig. 3 confirm this expectation. As can be seen from the displayed plots 
the general tendency for the changes caused by DM interaction is the same for XX and Heisenberg chains. 

To conclude, we have examined the stability of some spin-i transverse XX and Heisenberg chains with 
respect to dimerization in the presence of DM interaction analysing the dependence of the ground state energy 
on dimerization parameter. If the orientations of D„ are the same at all sites essential simplification occurs, 
i.e. the terms with cross product can be eliminated from the Hamiltonian resulting in the appearance of the 
anisotropy in the exchange interaction. In addition, the external field becomes more complicated having three 



site dependent components (see (Q)). The transformed Hamiltonian may be more transparent and convenient 
for further analytical or numerical treatment. For the transverse XX chain we have found that DM interaction 
having only z component may act both in favour of the dimerization or against it. The result of its influence 
depends on the dependence of DM interaction on the amplitude of lattice distortion in comparison with a 
corresponding dependence of the isotropic exchange interaction. For the Heisenberg chain DM interaction 
having the same orientation at all sites leads to XX anisotropy (and extra nonuniform on-site fields in the 
presence of an external field). For the Heisenberg chain without field we have observed qualitatively the same 
behaviour as in the XX chain determined by the dependence of DM interaction on the lattice distortion. 

O. D. and J. R. acknowledge the kind hospitality of the Max Planck Institute for the Physics of Complex 
Systems, Dresden in the end of 1999 when a main part of the paper was done. The present study was partly 
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FIGURE CAPTURES 

FIGURE 1. Dependence £(S) — £(0) @ vs S for the transverse XX chain with DM interaction. J = 2, 
f2o = 0, a = 0.8, D z = (curve 1), D z = 0.4, 0.8 for k = 1 (curves 2 and 3, respectively) and k = (curves 4 
and 5, respectively). Solid curves correspond to analytical calculations (N — > oo), whereas dotted ones to exact 
diagonalization results (N = 24). The dashed curve corresponds to the numerical results for N = 28. 



FIGURE 2. Dependence 6* vs a for the transverse XX chain with DM interaction. J = 2, a = 0.8, fio = 
(a), Oo = 0.2 (b), D z = (curve 1), Z? z = 0.4, 0.8 for k = 1 (curves 2 and 3, respectively) and k = (curves 4 
and 5, respectively). 

FIGURE 3. Dependence £{8) — £(0) vs S for the Heisenberg chain with DM interaction obtained by exact 
diagonalization of finite chains. J — 2, £1 = 0, a = 0.8, D = (curve 1), D = 0.4, 0.8 for k = 1 (curves 2 and 
3, respectively) and = (curves 4 and 5, respectively). We also plotted the dependence £{$) — £(0) vs 5 for 
the corresponding XX chain with D = (i.e. curves 1 in Fig. 1) as it follows from analytical calculation for 
N — > oo (dashed curve) and numerical computation for TV = 24 (dotted curve). 




Figure 1: Dependence £(S) — £ (0) (JSj) vs S for the transverse XX chain with DM interaction. J = 2, fio = 0, a = 0.8, 
_D Z = (curve 1), D z = 0.4, 0.8 for k = 1 (curves 2 and 3, respectively) and k = (curves 4 and 5, respectively). Solid 
curves correspond to analytical calculations (N — > oo), whereas dotted ones to exact diagonalization results (N = 24). 
The dashed curve corresponds to the numerical results for iV = 28. 




Figure 2: Dependence S* vs a for the transverse XX chain with DM interaction. J = 2, a = 0.8, fio = (a), fio = 
(b), D z = (curve 1), .D z = 0.4, 0.8 for k — 1 (curves 2 and 3, respectively) and k = (curves 4 and 5, respectively 
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Figure 3: Dependence £(S) —£(0) vs 8 for the Heisenberg chain with DM interaction obtained by exact diagonalization 
of finite chains. J = 2, flo = 0, a = 0.8, D — (curve 1), D = 0.4, 0.8 for k — 1 (curves 2 and 3, respectively) and 
k — (curves 4 and 5, respectively). We also plotted the dependence £(&) — £(0) vs 8 for the corresponding XX chain 
with D = (i.e. curves 1 in Fig. 1) as it follows from analytical calculation for iV — » oo (dashed curve) and numerical 
computation for N = 24 (dotted curve). 



